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Abstract 

Various nuclear structure observables are evaluated employing low-momentum nucleon- 
nucleon (NN) potentials Mow-k derived from the CD-Bonn and Nijmegen NN inter- 
actions Vnn. By construction, the high momentum modes of the original Vnn are 
integrated out in V\o^f,-k, with the requirement that the deuteron binding energy 
and low energy phase shifts of Vnn are exactly reproduced. Using this interaction, 
we evaluate the bulk properties (binding energy and saturation density) of nuclear 
matter and finite nuclei, in particular their dependence on the cut-off parameter. 
We also study the pairing gap and the residual interaction in nuclear matter in 
terms of the Landau parametrization. At low and medium densities, the HF and 
BHF binding energies for nuclear matter calculated with the l^ow-k's derived from 
the CD-Bonn and Nijmegen potentials are nearly identical. The pairing gaps and 
Landau parameters derived from Mow-k are remarkably close to those given by the 
full-space Vnn- The Viow-k interactions, however, fail to reproduce the saturation 
property of nuclear matter at higher densities if the cut-off for the high momentum 
modes is assumed density independent. 

PACS: 21.30.Fe; 13.75.Cs; 27.80.+j; 21.60.-n; 
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1 Introduction 



In the past decade, there has been much progress in treating low-energy nu- 
clear systems using the effective field theory (EFT) and renormalization group 
(RG) approaches, [1,2,3,4,5,6,7,8] and references therein. The fundamental 
theory of strong interaction, QCD, is non-perturbative in the energy regime 
characteristic for low energy nuclear physics. This poses a major difficulty for 
a proper derivation of the nuclear force from the first principles. However, this 
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difficulty can be overcome by the EFT-RG approach where one "integrates 
out" the fast fields beyond the chiral symmetry breaking scale ~ 1 GeV. 
Well below this scale, the appropriate degrees of freedom for nuclear systems 
are pions and nucleons. In this way, one obtains for these degrees of freedom 
a chiral symmetric EFT, which is used to construct either perturbative or 
non-perturbative controlled approximations for low energy nuclear systems. 

As is well known, there are several reahstic models for the NN potential Vatat, 
such as the CD-Bonn [9], Nijmegen [10] and the Argonne [11] potentials. These 
models all have the same long range tail, namely the one-pion exchange poten- 
tial. But the short range parts significantly differ from each other, despite the 
fact they all fit essentially exactly the empirical deuteron binding energy and 
low energy phase shifts up to Eiab ~350 MeV. Clearly the short range (high 
momentum) parts of these potentials are model dependent and are different 
e.g. for the local and non-local interaction models. Employing the EFT-RG 
ideas, Refs. [12,13,14] "integrated out" the high momentum modes of these 
potentials and thereby obtained a low momentum NN potential Viow-k- The 
high momentum modes of the original Vjvat are eliminated in Viow-k with the 
requirement that the deuteron binding energy and low energy phase shifts of 
Vnn are exactly reproduced. This is in line with the central theme of the EFT- 
RG approach that physics in the infrared region is insensitive to the details of 
the short distance dynamics. 

The purpose of this work is to calculate various nuclear structure observables 
for infinite nuclear matter and finite nuclei starting from the Viow-k interaction 
derived for the CD-Bonn and Nijmegen potentials [12,13,14]. In particular, we 
would like to explore to which extent the results derived from l^ow-k differ 
from the predictions obtained for the underlying interaction model Vnn and 
the sensitivity of the results on the cut-off parameter A. For that purpose we 
evaluate the ground-state energy of nuclear matter and finite nuclei taking 
^^O as a typical example. Apart from the bulk properties of nuclear systems 
we also investigate the predictions for the residual interaction characterized 
by the pairing gap A of neutron-neutron or proton-proton pairing as well 
as the Landau parametrization of the particle-hole interaction for low-energy 
excitations in nuclear matter. 

The paper is organized as follows. In section 2 we shortly review the con- 
struction of the low-momentum Viow-k ■ The methods for the evaluation of 
ground-state properties of nuclear matter and finite nuclei as well as the cor- 
responding results are presented in sections 3 and 4, respectively. Section 5 
contains a discussion of the gap equation and the Landau parameters charac- 
terizing the residual interaction. The conclusions of this study are summarized 
in section 6. 
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2 Determination of low momentum nucleon-nucleon potentials 



We start with a brief review of the method of carrying out the high- momentum 
decimation [12,13,14]. In accordance with the general definition of a renormal- 
ization group transformation, the decimation must be such that low energy 
observables calculated from the full Vnn must be reproduced by V|ow-kj but 
with all loop integrals cut off at some A. To meet this requirement, a T-matrix 
equivalence approach is used [12,13,14]. One starts from the half-on-shell T- 
matrix 

T{k',k,k') = VMk',k)+ q'dqVMk', q) ^2 _ ^2 + iQ+ T{q, k, k'), (1) 
and then defines an effective low-momentum T-matrix by 

71ow-k(p',P,P^) = Mow-k(p',p)+ q^dqViov,--k{p' , g)^^-^— ^riow-k(?,P,P^), 

(2) 

where A denotes the momentum space cut-off. The above T-matrices satisfy 
the condition 

7^(p',P,p') = Tiow-k(p',P,p'); for (p',p)<A. (3) 

Viow-k has been obtained by solving the above equations using the procedures 
of Ref. [12,13]. In particular the effects of the folded diagrams are taken into 
account to derive an effective interaction Viow-k which is energy-independent. 

We have performed numerical checks to make sure that the resulting Viow-k 
does preserve the deuteron binding energy and low energy phase shifts of 
the original NN potential. This has been repeated assuming various cut off 
parameters A using the CD Bonn [9] and the Nijmegen [10] interaction models 
for the bare interaction Vnn- For both of these interactions we have ignored the 
charge dependence and used only the proton-neutron part for all interactions 
channels. 



3 Energy of Nuclear Matter 



It is well known that the strong short-range and tensor components of re- 
alistic NN interactions induce strong NN correlations in the wave function 
of nuclear many-body system. Therefore rather sophisticated approximation 
schemes for the solution of the many-body problem are required to obtain rea- 
sonable results for nuclear structure observables. Such approximation schemes 
must account at least for the two-nucleon correlations. Typical examples are 
variational approaches using correlated basis functions [15], the coupled cluster 



3 



approach [16], the self-consistent evaluation of Green's functions [17] and the 
hole-line expansion within the Brueckner-Hartree-Fock (BHF) approximation 
[18]. 

The effective interaction Viow-k has been constructed to account for NN cor- 
relations with momenta larger than the cut-off momentum A as they occur 
for two interacting nucleons in the vacuum. The question which we now want 
to address is whether such a treatment of correlations is sufficient to allow 
for simple mean field or Hartree-Fock calculations of the nuclear many-body 
system. As a first example we consider the equation of state of infinite nuclear 
matter. 

The left panel of Fig. 1 shows the results of Hartree-Fock (HF) calculations for 
the energy per nucleon of symmetric nuclear matter as a function of density, 
which is characterized in this figure by the corresponding Fermi momentum 
kp- For the bare CD-Bonn interaction [9] the HF calculations yield positive 
energies ranging between 3 MeV per nucleon at small densities and 20 MeV 
per nucleon for kp around 1.8 fm^^ which corresponds to a density of around 
2.3 times the empirical saturation density. This demonstrates the well known 
fact that even modern non-local models for a realistic NN interaction, which 
tend to be "soft" , yield unbound nuclei if the effects of correlations are ignored 
[19]. 

This lack of binding energy can be cured if the effects of NN correlations are 
taken into account, e.g. by means of the BHF approximation. In the BHF 
approximation one replaces the bare interaction by the G-matrix which obeys 
the Bethe-Goldstone equation. For a pair of nucleons with vanishing center of 
mass momentum it takes the form 

G{k',k,2ek)=VNN{k',k)+ / q^dqVNN{k',q)- ^—^G{q,k,26k). 

Jo ISk — ^Eq -\- lU^ 

(4) 

Comparing this equation with the Lippmann-Schwinger equation (1), one can 
see that the G-matrix is closely connected to the T-matrix except for the 
fact that the Pauli operator Q in (4) prevents scattering into states which 
are occupied by other nucleons and that the propagator for the two nucleons, 
which for the vacuum is defined in terms of the kinetic energies (see eq. (1)) is 
replaced in (4) by the propagator in the medium defined in terms of the single- 
particle energies e^. Both of these modifications tend to reduce the effects of 
correlations and are referred to as Pauli and dispersive quenching, respectively 
[20]. 

It should be mentioned that we employ a very simple version of the BHF ap- 
proach, in which the single-particle energies for states below the Fermi energy 
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are parametrized in terms of an effective mass and a constant energy shift 



while for states above the Fermi energy the single-particle energies are replaced 
by the kinetic energy (the so-called conventional choice) . The Pauli operator is 
treated in the angle-averaged approximation. A more sophisticated treatment 
of the two-nucleon propagator in the Bethe-Goldstone equation leads to larger 
binding energies [21]. Note, however, that it is not the aim of this study to 
optimize the BHF approximation, but to study the possible use of the Viow-k 
interaction in restricted model-space calculations. 

Results of such BHF calculations are displayed in the right panel of Fig. 1. 
The conventional BHF calculation for the bare CD Bonn interaction yields 
an energy of -17 MeV per nucleon and a saturation density around twice 
the empirical saturation density. The left panel of Fig. 1 also shows the HF 
energies calculated for Viow-k assuming three different values for the cut-off 
parameter A. The resulting energies are attractive. This demonstrates that the 
correlations from two-particle states with momenta q larger than the cut-off 
A, which are accounted for in Mow-k, are sufficient to provide this attraction. 
From this consideration it is evident that the energies are more attractive the 
smaller the cut-off A, a feature seen in Fig. 1. 

Next, in performing BHF calculations for Viow-k, we must take into account 
that V[ow-k is designed for a model space with relative momenta smaller than 
A, i.e. the matrix elements of VJow-k must be supplemented by the value zero 
for momenta q larger than the cut-off. In other words, the integral in the 
Bcthe-Goldstone eq. (4) is restricted to momenta q below the cut-off in order 
to avoid double counting of the contributions from momenta higher than the 
cut-off. 

The energies resulting from such BHF calculations (see left panel of Fig. 1) 
arc below the corresponding HF energies as additional correlations from NN 
states above the Fermi momentum but below the cut-off are taken into ac- 
count. One also observes that the BHF calculations using VJow-k lead to more 
attractive energies than the BHF calculation using the bare interaction. This 
is due to the fact, that the dispersive quenching effects in the G-matrix dis- 
cussed above are taken into account only for intermediate two-particle states 
with momenta below the cut-off A, while the BHF calculation using the bare 
interaction considers these effects for all states. Since the dispersive quenching 
is a very important ingredient of BHF calculations, that provides the satura- 
tion effect, the BHF calculations employing Vjow-k do not exhibit a minimum 
in the energy versus density curve. One may conclude that Viow-k yields a good 
approximation for BHF calculations at small densities. In this regime even a 
Hartree-Fock calculation leads to reasonable results, if A is small enough. At 
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higher densities, however, Viow-k leads to too much binding as it does not 
account for the dispersive quenching of high-momentum modes. 

In Fig. 2 we compare results derived from two different NN interaction mod- 
els: the CD-Bonn (used above) and the Nijmegen [10] potentials. One can see 
from this figure that the HF calculations using the Nijmegen interaction yield 
energies which are more repulsive than those derived from the CD-Bonn inter- 
action. The non-local CD-Bonn is "softer" than the local Nijmegen interaction. 
As it has been discussed e.g. in Ref. [19] the softer interaction model yields 
more attractive energies also within the BHF approximation (see right panel 
of Fig. 2). Although the results of HF as well as BHF calculations arc rather 
sensitive to the choice of the interaction, the corresponding results from Viow-k 
arc rather insensitive to the interaction model on which they arc based. This 
is true for both the HF and the BHF approximations. This can be understood 
from the fact that the T-matrix is essentially independent on the interaction 
used, as all realistic interactions arc fitted to the empirical phase shifts. Since 
the definition of Vlow-k is related to the T-matrix (see eq. (3)), the resulting 
effective interaction is rather insensitive to the original interaction it is based 
on. 



4 Finite Nuclei 



Since the use of Viow-k yields rather reasonable results for the energy of nuclear 
matter at small densities it is of interest to investigate the predictions of Viow-k 
for the ground-state properties of light nuclei. For that purpose we employ the 
calculation scheme, which has recently been discussed in ref. [22] and apply 
it for the nucleus ^^O. This scheme is an extension of the conventional BHF 
approximation for finite nuclei. It assumes a model space, which in our example 
is defined in terms of configurations of nucleons occupying oscillator states 
up to the IpOf shell. Correlations outside this model space, which we would 
like to call short-range correlations, are treated within the framework of the 
BHF approximation. The long-range correlations from configurations inside 
the model space are treated within the framework of the Green's function 
approach, assuming a self-energy which accounts for 2-particle - 1-hole and 
2-hole - 1-particle admixtures. 

This Green's function approach (for further details see ref. [22]) yields the 
spectral function of the hole-state distribution S'^ij{uj), which is the probability 
to remove a nucleon with quantum numbers r for isospin, / for orbital angular 
momentum and j for the total angular momentum by absorption of an energy 
—CO. From this spectral function we can easily determine the total occupation 
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probability 



(6) 



the mean single-particle energy of the distribution function 



J —oo 



(7) 



and the total energy 



Tlj 



(2j + 1) 



2 



— oo 




) 
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Results for these observablcs are listed in table 1 using the bare CD-Bonn 
interaction along with Viow-k derived from the CD-Bonn potential with A = 2 
fm~^ and 3 fm~^. For the interactions Mow-k we have performed calculations 
in which the effects of short-range correlations are ignored (labeled HF) and 
taken into account (labeled BHF). 

The BHF calculations employing V^ow-k yield results which are very close to 
the corresponding ones for the bare interaction. The agreement is better for 
the larger value of the cut-off A. The BHF calculations using V\ow-k yield 
energies which are more attractive than the corresponding result for the bare 
interaction. This can be explained in the same way as for the case of nuclear 
matter: V[ow-k accounts for correlations originating from states above the cut- 
off as they occur for two free nucleons. The dispersive quenching effect in the 
G-matrix calculated for the bare interaction, however, reduces these attractive 
correlation contributions. We find that this effect is rather small for finite 
nuclei, comparable to the same effect in nuclear matter at small densities (see 



The predictions from calculations using Viow-k in the HF approach yield en- 
ergies which are less attractive than those from the bare interaction. This is 
due to the fact that they do not account for short-range correlations originat- 
ing from configurations outside the model space (up to IpOf shell) but only 
below the cut-off. Nevertheless, the results obtained for A = 2 fm"-*^ are quite 
reasonable. 

The predictions for the occupation probabilities are larger for the calculations 
using V^iow-k than for the bare interaction V. This is due to the fact that the HF 
calculations using Mow-k only account for a depletion of the occupation due to 
long-range correlations within the IpOf model space, while the corresponding 
BHF calculations include effects of correlations only up to the cut-off A. One 
may observe, however, that the BHF calculation using Viow-k with A = 3 fm~^ 
yields results for the occupation which are only slightly above those derived 
from the bare interaction. 



Fig. 1). 
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The reliability of Viow-k in low-energy nuclear structure calculations of finite 
nuclei is also supported by the results for the spectral function displayed in 
Fig. 3. The structure of the spectral function derived from the bare interaction, 
represented by the dashed line in both parts of the figure, is nicely reproduced 
by the spectral function derived from Viow-k- This is true in particular for 
the BHF approach, displayed by the solid line in the upper panel, but also 
holds for the HP approximation, which is shown in the lower panel. The main 
deficiency is a shift of the strength to larger energies a;, a feature which we 
have already outlined above in discussing the mean single-particle energies 
shown in table 1. 



5 Pairing and Residual Interaction in Nuclear Matter 

Up to this point we have focused the discussion on observables, which are 
characteristic for the ground state of nuclear systems. In this section we would 
like to address the question how reliable is the use of Viow-k for observables 
related to the low-energy excitation spectrum. As the first example we consider 
the evaluation of the pairing gap in nuclear matter. In particular we take the 
example of neutron- neutron pairing in the ^5*0 channel for neutron matter at 
kp — 1.36 fm""^. This pairing gap is determined from the solution of the gap 
equation (see e.g. [23,24]) 



Here we will again parametrize the single-particle spectrum by eq. (5) and 
consider the cases m* — m and m* — 700 MeV. For a realistic interaction, 
one has to integrate up to rather large values of k' before a convergent result 
is obtained. Corresponding results for A{k) for the bare CD-Bonn interaction 
are displayed in Fig. 4. For the Viow-k interactions, the integral in (9) can be 
constrained to values of k' below the cut-off A. 

The agreement for the resulting gap functions is very good for momenta below 
the cut-off A and a single-particle spectrum with an effective mass identical to 
the bare mass. In this case the use of Viow-k in the gap equation corresponds to 
the use of the effective model-space interaction of ref. [24]. However, if we use 
a single-particle spectrum renormalized by the medium, the results obtained 
from V\ovj-k differ in a more significant way from the corresponding result 
obtained for the bare CD-Bonn interaction. If, e.g. we use an effective mass 
m* of 700 MeV (see right panel of Fig. 4), Viow-k tends to overestimate the 
value for the gap A as it does not account for the reduction of the contributions 
from high lying NN states, which is due to the decrease of the effective mass. 



2 r 

TT Jo 



dk' k'^ V{k, k') 



A{k') 



(9) 



2v/(4 - ^F? + A{k'y 



8 



This feature is identical to the effect of the dispersive quenching of correlations 
discussed above. 



Finally, we would like to consider the residual interaction for low-energy ex- 
citations as characterized by the Landau parameters. It is not our goal to 
employ a very sophisticated calculational scheme for the Landau parameters. 
This would involve in particular a careful analysis of the contributions due 
to the induced interaction terms [25,26,27]. Our aim here is to compare the 
results from the bare interaction Vnn, the corresponding G-matrix and those 
derived from l^ow-k, in order to see to which extent Mow-k accounts for the 
short-range correlation effects in a nuclear system. For that purpose, we deter- 
mine the Landau parameters either from the bare interaction or the G-matrix 



The results of our calculations of the Landau parameters are listed in table 2. 
Beside those parameters we also list the effective mass determined by 

^ = i + k (10) 

m 3 

and the coefficient /3 for the symmetry energy 

Comparing the various approximations with our "best" result derived from 
the G-matrix of the CD-Bonn potential, one finds that the bare interaction 
Viow-k with a small cut-off parameter yields a good approximation if one would 
like to restrict oneself to the bare interaction. If an iteration is carried out in 
the particle-particle channels for momenta up to the cut-off A, one obtains the 
best results using a larger cut-off. 



6 Conclusions 



The aim of this study was to explore to which extent the effects of short- 
range correlations, which are due to excited states above a cut-off parameter 
A, are accounted for by the effective interaction Viow-k? which has recently 
been constructed [12,13,14]. This effective interaction accounts for short-range 
correlations of two-nucleons interacting in the vacuum; the question is whether 
this is a good approximation also for the corresponding correlations of nucleons 
interacting in a nuclear medium. To study this question we considered various 
nuclear structure observables. 

Inspecting the equation of state calculated for nuclear matter one may con- 
clude that Mow-k yields a good approximation for BHF calculations at small 
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densities. In this regime even a Hartree-Fock calculation leads to reasonable 
results, if A is small enough. At higher densities, however, Viow-k produces 
too much binding as it does not account for the dispersive quenching effects, 
which are taken into account e.g. in the G-matrix. Therefore calculations using 
Viow-k cannot reproduce the saturation feature of nuclear matter. 

The Viow-k interaction is well suited for studies of finite nuclei like ^^O, which is 
again a system of small density. In this case we obtain reasonable agreement for 
the calculated energies, but also for the spectral function for nucleon removal. 

Furthermore, we demonstrated that the use of Mow-k is also very efficient in 
calculating quantities like the pairing gap in infinite nuclear matter. The agree- 
ment with the exact result is very good if one may replace the single-particle 
energy above the cut-off by pure kinetic energy. Deviations are observed if this 
is not the case. These deviations are again related to the dispersive quenching 
mechanism mentioned above. 

In conclusion we would like to emphasize that Viow-k can be a very useful 
tool for low-energy nuclear structure calculations. However, one should be 
cautious if the observable of interest is sensitive to the single-particle spectrum 
at energies above the cut-off. This is a problem in particular for high densities. 
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HF 






BHF 








A=2 Im-i A= 


=3 fm-i 


A=2 fm-i 


A=3 fm-i 


bMro V 


Exp 


£sl/2 


-39.28 


-31.31 


-45.23 


-44.74 


-44.00 


-40 ± 8 


ep3/2 


-18.32 


-12.34 


-25.07 


-24.83 


-24.29 


-18.45 




-15.01 


-9.01 


-21.97 


-21.81 


-21.26 


-12.13 


Psl/2 


0.964 


0.964 


0.906 


0.880 


0.867 




Pp3/2 


0.928 


0.936 


0.844 


0.828 


0.823 




Ppl/2 


0.912 


0.912 


0.832 


0.815 


0.802 




E/A 


-4.96 


-1.69 


-8.22 


-8.08 


-7.78 


-7.98 



Table 1 

Single-particle energies (e^) and occupation probabilities (pi) for protons in ^^O and 
the total energy per nucleon (E/A) as calculated from various approximations and 
interactions are compared to the experimental data. The CD Bonn interaction has 
been used (bare V) and the effective interactions V|ow-k derived from it assuming 
a cut-off A of 2.0 fm~^ and 3.0 fm~^. The effects of short-range correlations have 
been ignored (HF) or taken into account by means of the BHF approximation. All 
energies are given in MeV. 







V 






G 






A=2 fm-i 


A=3 fm-i 


bare V 


A=2 fm-i 


A=3 fm-i 


bare V 


/o 


-1.819 


-1.315 


-0.465 


-2.003 


-1.807 


-1.601 




0.502 


0.235 


0.072 


0.614 


0.492 


0.422 


90 


0.250 


0.362 


0.226 


0.211 


0.279 


0.276 




0.970 


0.899 


0.709 


0.998 


0.973 


0.929 


ho 


1.788 


1.759 


1.618 


1.753 


1.717 


1.694 


K 


-0.596 


-0.631 


-.581 


-.588 


-.620 


-.613 


m* /m 


0.702 


0.684 


0.620 


0.696 


0.688 


0.676 


P [MeV] 


28.98 


24.46 


23.43 


31.40 


29.35 


28.48 



Table 2 



Landau parameters, effective mass and the coefficient for the symmetry energy f3 as 
calculated from the interaction V and the corresponding G matrix in nuclear matter 
at a Fermi momentum kp of 1.36 fm~^. The CD Bonn interaction has been used 
(bare V) and the effective interactions VJow-k derived from it assuming a cut-off A 
of 2.0 fm^^ and 3.0 fm^^. The Landau parameters are normalized with respect to 
the density of states at the Fermi momentum using an effective mass as listed. 
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Fermi Momentum [fm ] Fermi Momentum [fm ] 

Fig. 1. Binding energy of nuclear matter as a function of the Fermi momentum for 
tire HF (left panel) and BHF approximations (right panel). The bare interaction 
(CD-Bonn [9]) and corresponding Viow-k interactions with various cut-off parame- 
ters A were used. 
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Fig. 2. Binding energy of nuclear matter as a function of the Fermi momentum 
for the HF (left panel) and BHF approximations (right panel). The interactions 
are modeled by the CD-Bonn and Nijmegen potentials and corresponding Viow-k 
interactions with a cut-off parameter A = 2.5 fm~^. 
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Fig. 3. The spectral function for proton hole strength in ^^O in the channel 
for the BHF (upper panel) and the HF (lower panel) approximations. The spectral 
distribution shown by the solid line has been calculated for Vjow-k with A = 2 fm~^. 
For comparison, we show the results derived from the BHF approximation using the 
bare interaction (dashed lines) [22]. 
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Fig. 4. The pairing gap A(A;) for neutron-neutron pairing in the ^Sq channel for 
neutron matter with a Fermi momentum kp = 1.36 fm~^ as function of the nucleon 
momentum k. The results obtained for the bare CD-Bonn interaction (solid lines) are 
compared to those derived from Viow-k using cut-offs A = 2 fm~^ and 3 fm^^. The 
left and right panels show the results for non-renormalized single-particle spectrum 
(m* = m) and for the renormalized spectrum with an effective mass m* = 700 MeV, 
respectively. 
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